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[i]  Global  comparisons  of  barotropic  and  internal  tides  generated  in  an  eddy-resolving 
ocean  circulation  model  are  made  with  tidal  estimates  obtained  from  altimetric  sea  surface 
heights  and  an  altimetry-constrained  tide  model.  As  far  as  we  know,  our  Hybrid  Coordinate 
Ocean  Model  (HYCOM)  simulations  shown  here  and  in  an  earlier  paper  are  the  only 
published  high-resolution  global  simulations  to  contain  barotropic  tides,  internal  tides,  the 
general  circulation,  and  mesoscale  eddies  concurrently.  Comparing  the  model  barotropic 
tide  with  a  global  data-assimilative  shallow  water  tide  model  shows  that  the  global  tidal 
elevation  differences  are  approximately  evenly  split  between  discrepancies  in  tidal 
amplitude  and  phase.  Both  the  model  and  observations  show  strong  generation  of  internal 
tides  at  a  limited  number  of  “hot  spot”  regions  with  propagation  of  beams  of  energy  for 
thousands  of  kilometers  away  from  the  sources.  The  model  internal  tidal  amplitudes 
compare  well  with  observations  near  these  energetic  tidal  regions.  Averaged  over  these 
regions,  the  model  and  observation  internal  tide  amplitude  estimates  agree  to 
approximately  15%  for  the  four  largest  semidiurnal  constituents  and  23%  for  the  four 
largest  diurnal  constituents.  Away  from  the  hot  spots,  the  comparison  between  the  model 
and  altimetric  amplitude  is  not  as  good  due,  in  part,  to  two  problems,  errors  in  the  model 
barotropic  tides  and  overestimation  of  the  altimetric  tides  in  regions  of  strong  mesoscale 
eddy  activity.  Examining  the  general  energy  distribution  of  the  simulated  internal  tide  is  an 
important  first  step  in  the  evaluation  of  internal  tides  in  HYCOM. 

Citation:  Shriver,  J.  F.,  B.  K.  Arbic,  J.  G.  Richman,  R.  D„  Ray,  E.  J.  Metzger,  A.  J.  Wallcraft,  and  P.  G.  Timko  (2012), 

An  evaluation  of  the  barotropic  and  internal  tides  in  a  high-resolution  global  ocean  circulation  model,  J.  Geophys.  Res.y  117, 
C 1 0024,  doi:  10.1 029/20 1 2 JC008 1 70. 


1.  Introduction 

[2]  Oceanic  internal  tides  are  internal  waves  with  tidal 
periodicity  that  are  generated  by  the  interaction  of  barotropic 
tidal  currents  with  variable  bottom  topography.  Freely 
propagating  diurnal  internal  tides  are  theoretically  con¬ 
strained  to  the  approximate  latitude  range  30°S-30°N  [Gill, 
1982,  p.  258],  with  semidiurnal  tides  found  in  this  range 
and  higher  latitudes.  They  play  a  key  role  in  dissipating  tidal 
energy  and  mixing  in  the  deep  ocean  [e.g.,  Egbert  and  Ray , 
2000;  Jayne  and  St.  Laurent ,  2001 ;  Simmons  et  al. ,  2004]. 

[3]  A  first  attempt  to  resolve  internal  tides,  along  with 
barotropic  tides  and  the  eddying  general  circulation,  in  a 
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global  ocean  circulation  model  is  described  by  Arbic  et  al. 
[2010].  The  simulations  in  the  work  by  Arbic  et  al.  [2010] 
and  in  this  study  utilize  the  Hybrid  Coordinate  Ocean 
Model  (HYCOM),  which  is  being  developed  by  the  United 
States  Navy  as  an  ocean  nowcast/forecast  model.  These 
simulations  employ  both  tidal  and  atmospheric  forcing, 
in  contrast  to  earlier  global  baroclinic  tide  simulations  that 
employed  only  tidal  forcing  [Arbic  et  al ,  2004;  Simmons 
et  al .,  2004].  Atmospheric  forcing  allows  for  a  more  realis¬ 
tic  horizontally  varying  stratification  in  HYCOM,  in  contrast 
to  the  uniform  stratification  employed  in  the  work  by  Arbic 
et  al.  [2004]  and  Simmons  et  al.  [2004] 

[4]  The  HYCOM  tides  simulations,  which  thus  far  have 
been  run  only  in  forward  (nonassimilative)  mode,  are  contin¬ 
ually  being  updated,  especially  with  regards  to  the  parame¬ 
terized  topographic  wave  drag.  As  a  result,  the  simulations 
utibzed  here  are  not  identical  to  those  by  Arbic  et  al  [2010], 
though  the  accuracy  of  the  barotropic  tide  is  comparable.  In  the 
simulation  analyzed  here  we  capture  93.2%  of  the  sea  surface 
elevation  variance  of  the  eight  largest  tidal  constituents  in  the 
standard  set  of  102  pelagic  tide  gauges  [Shum  et  al.,  1997]; 
in  the  work  by  Arbic  et  al.  [2010]  they  captured  92.6%. 

[5]  As  far  as  we  know,  our  HYCOM  simulations  are  the 
only  published  global  simulations  to  contain  barotropic 
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Figure  1.  Standard  errors  (cm)  of  altimeter-based  along-track  estimates  of  the  Oj  surface  tidal  elevation 
(a)  without  and  (b)  with  prior  correction  of  the  altimetry  for  nontidal  sea  surface  variability.  The  largest 
errors  in  Figure  la  reach  3  cm.  Results  for  M2  are  very  similar.  Approximately  17  years  of  T/P  and  Jason 
altimetry  are  used  in  the  tidal  estimation.  See  also  Carrere  et  al  [2004,  Figure  5]. 


tides,  internal  tides,  the  general  circulation,  and  mesoscale 
eddies  concurrently  and  at  high  horizontal  resolution.  This 
paper  examines  how  accurately  HYCOM,  forced  only  by 
atmospheric  forcing  and  the  astronomical  tidal  potential,  can 
simulate  the  global  barotropic  tide  fields  in  comparison  to  an 


altimetry-constrained  barotropic  tide  model  [Egbert  et  al. , 
1994]  and  the  global  internal  tide  fields  in  comparison  to 
an  observed  data  set.  The  only  global  observations  of  baro¬ 
tropic  and  internal  tides  are  based  on  satellite  altimetiy.  We 
compare  the  barotropic  tides  in  HYCOM  to  output  from  an 
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Figure  2.  Amplitude  (cm)  of  M2  surface  tidal  elevation  in  (a)  TPX07.2  (an  update  to  that  described  by 
Egbert  et  al  [1994]),  a  barotropic  tide  model  constrained  by  satellite  altimetiy,  and  (b)  HYCOM  simula¬ 
tions  in  which  the  tide  is  unconstrained  by  satellite  altimetry.  Lines  of  constant  phase  plotted  every  45°  in 
Figures  2a  and  2b  are  overlaid  in  white. 


altimetiy-constrained  barotropic  tide  model  [Egbert  et  aly 
1994]  and  the  internal  tides  in  HYCOM  to  results  from  an 
analysis  of  along-track  satellite  altimetry  data  [Ray  and 
Mitchum ,  1996].  Several  previous  comparisons  of  modeled 


and  observed  internal  tides  have  utilized  regional  models  of 
strong  internal  tide  generation  sites  forced  by  specified  bar¬ 
otropic  tides  at  their  horizontal  boundaries  [e.g.,  Cummins 
et  al.y  2001;  Kang  et  al.y  2000;  Merrifield  et  al.y  2001]. 
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Figure  3.  Same  as  Figure  2  but  for  the  Kj  tidal  constituent. 


Arbic  et  al  [2010]  validated  the  HYCOM  barotropic  tides 
via  comparison  to  a  standard  set  of  pelagic  tide  gauges 
[Shum  et  al,  1997],  which  obviously  offer  much  less  uni¬ 
form  coverage  of  the  global  ocean  than  satellite  altimetry. 
The  validation  of  internal  tides  by  Arbic  et  al  [2010]  was 
done  for  only  one  tidal  constituent  (M2)  and  for  a  limited 
area  around  Hawaii.  Our  results  represent  the  first  global 


quantitative  comparison  of  the  simulated  internal  tide  field  to 
an  observed  data  set. 

2.  Model  and  Data 

[6]  As  in  the  work  by  Arbic  et  al  [2010],  the  HYCOM 
simulation  examined  in  this  study  utilizes  gcopotential  tidal 


C10024 


SHRIVER  ET  AL.:  BAROTROPIC  AND  INTERNAL  TIDES  IN  HYCOM 


C10024 


Table  1.  Global-Averaged  Amplitude,  Phase,  and  Total  RMS 
Errors  of  HYCOM  Surface  Tidal  Elevations  Measured  Against 
TPX07.2,  Following  the  Error  Derivation  in  (3)3 


Amplitude  Error 

Phase  Error 

Total  Error 

m2 

5.560 

4.999 

7.477 

S2 

2.807 

3.544 

4.522 

N2 

0.810 

0.794 

1.134 

k2 

0.770 

1.166 

1.397 

K, 

1.555 

1.627 

2.251 

0, 

1.714 

1.467 

2.257 

Pi 

0  498 

0.543 

0.737 

Qi 

0.379 

0.371 

0.530 

*Units  are  in  centimeters. 


foremg  for  M2,  S2,  N2,  and  K2  (the  four  largest  semidiurnal 
constituents),  and  Oj,  Pi,  and  Qi  (the  four  laigest  diurnal 
constituents),  a  scalar  self-attraction  and  loading  correction 
(SAL)  [Ray,  1998],  a  parameterized  topographic  wave  drag, 
32  layers  in  the  vertical  direction  and  a  nominal  horizontal 
resolution  of  1/12.5°  at  the  equator.  For  the  tides,  there  are 
only  two  adjustable  parameters,  the  scalar  SAL  and  the 
topographic  drag  amplification  factor.  These  two  parameters 
are  adjusted,  using  a  one-layer  barotropie  M2  only  version  of 
the  model,  to  minimize  the  differences  between  the  model 
M2  tide  and  the  102  pelagic  tide  gauges.  All  other  para¬ 
meters  are  the  same  as  the  parameters  in  the  nontidal  global 
model  at  the  same  resolution.  Therefore,  no  further  tuning 
was  done  specifically  for  the  internal  tides  in  the  model. 
Arbic  et  al  [2010]  describe  the  necessity  for  a  parameterized 
topographic  wave  drag  in  global  baroclinie  tide  models.  In 
the  simulation  by  Arbic  et  al  [2010],  the  topographic  wave 
drag  from  Arbic  et  al  [2004]  is  multiplied  by  a  factor  of  6 
giving  e-folding  time  scales  from  the  drag  of  1.5  h  to  6  days 
with  no  topographic  wave  drag  over  73%  of  the  world 
ocean.  To  mitigate  instability  from  the  extremely  short  e- 
folding  time  scales,  the  wave  drag  e-folding  time  in  our 
simulation  is  clipped  at  24  h  and  the  sealing  factor  increased 
to  1 2  to  keep  the  average  wave  drag  the  same  as  in  the  work 
by  Arbic  et  al  [2010].  For  additional  details  on  global  eddy- 
resol  vmg  HYCOM  the  reader  is  referred  to  Metzger  et  al 
[2010]. 

[7]  The  model  was  run  mterannually  over  the  period  7/2003- 
12/2010  using  3-hourly  Fleet  Numerical  Meteorology  and 
Oceanography  Center  Navy  Operational  Global  Atmospheric 
Prediction  System  (FNMOC  NOGAPS)  [Rosmond  et  al , 
2002]  atmospheric  forcing  with  wind  speeds  scaled  to  be  con¬ 
sistent  with  QuikSCAT  observations.  Total  sea  surface  height 
(SSH)  snapshots  were  saved  once  per  hour  for  this  period.  Since 
the  majority  of  low  vertical  mode  internal  tide  energy  in 
altimetry  is  thought  to  be  stationary  [Ray  and  Zaron ,  2011], 
we  used  model  results  from  calendar  year  2006  for  this  study. 

[8]  The  HYCOM  tidal  sea  surface  elevation  amplitude  and 
phases  were  calculated  as  a  complex  amplitude  using  standard 
harmonic  analysis  [Foreman,  1977]  applied  to  the  HYCOM 
total  SSH.  The  HYCOM  tidal  sea  surface  elevations  are 
dominated  by  the  barotropie  tides,  and  are  compared  to  a 
hydrodynamic  model  of  the  barotropie  tides  constrained  by 
satellite  altimetry  (TPX07.2,  an  update  to  that  described  by 
Egbert  et  al  [1994]).  In  this  comparison,  referred  to  hereafter 
as  the  “barotropie”  comparison,  the  HYCOM  results  are 
interpolated  to  the  lower-resolution  TPXO  grid. 


[9]  The  internal  tides  in  HYCOM  are  compared  to  alti- 
metric  internal  tidal  estimates  derived  from  approximately 
17  years  of  along-traek  TOPEX7POSEIDON  and  Jason  sat¬ 
ellite  altimetry.  To  facilitate  the  comparison  of  internal  tides 
in  HYCOM  with  those  in  along-track  altimetry  data  (here¬ 
after,  the  “internal  tide”  comparisons),  the  complex  ampli¬ 
tudes  from  the  HYCOM  tidal  analysis  are  first  interpolated  to 
the  along-track  altimeter  data  locations.  The  internal  tide 
complex  amplitudes  are  recovered  from  the  HYCOM  and 
along-traek  altimeter  analyses  via  band-pass  filtering  to  per¬ 
mit  wavelengths  in  the  50-400  km  range.  This  range  spans 
the  length  scales  of  the  low-mode  internal  waves  that 
HYCOM  is  able  to  resolve.  All  filtering  and  interpolating  is 
done  in  complex  spaee,  with  amplitude  discussed  in  seetion  3 
representing  the  positive  definite  magnitude.  All  analyses 
discussed  in  this  paper  foeus  on  locations  where  the  seafloor 
depth  exceeds  1500  m. 

[10]  For  the  satellite  data,  a  response  analysis  [Cartwright 
and  Ray ,  1990]  is  used  for  the  diurnal  and  semidiurnal 
bands,  supplemented  with  estimates  of  the  annual  cycle  and  a 
single  quarter-diurnal  constituent.  Solid  earth  tides  (includ¬ 
ing  the  component  arising  from  crustal  loading)  are  removed 
via  models.  Tides  are  independently  estimated  point-by-point 
along  the  satellite  tracks  [Ray  and  Mitchum ,  1996]. 

[11]  As  is  well  known,  the  satellite  repeat  period 
(9.9156  days)  aliases  all  diurnal  and  subdiumal  tides  to  long 
periods:  roughly  60  days  for  the  two  largest  semidiurnal 
tides  and  1 73  days  for  Kj .  The  alias  periods  for  the  tides  are 
similar  to  the  mesoseale  eddy  variability  time  scales,  and  the 
spatial  scales  for  low  vertical  mode  internal  waves  of  tidal 
period  are  similar  to  the  spatial  scales  of  mesoseale  eddies. 
Real  ocean  variability  at  these  alias  periods  can  direetly 
corrupt  tidal  estimation  (“mesoseale  contamination”)  when 
attempted  at  single  locations  along  track. 

[12]  In  fact,  extraction  of  internal  tide  signals  from  altimetry 
is  especially  problematic  in  regions  of  high  mesoseale  vari¬ 
ability  [Tierney  et  al ,  1998;  Carrere  et  al ,  2004].  This  diffi¬ 
culty  is  reflected  in  our  estimation  standard  errors,  shown  in 
Figure  la  for  the  Oi  constituent.  It  is  clear  that  the  largest 
errors  are  associated  with  strong  boundary  currents.  Although 
the  standard  errors  shown  here  formally  pertain  to  all  wave 
numbers,  the  spatial  scales  of  boundary  currents  can  drive 
these  largest  errors  into  the  internal  tide  band.  The  magnitudes 
of  the  errors  can  greatly  exceed  expected  internal  tide  ampli¬ 
tudes  (see  section  3.2),  making  signal  extraction  in  these 
regions  very  difficult. 

[13]  The  problem  of  mesoseale  contamination  can  be 
reduced,  but  not  eliminated,  by  removing  from  the  altimetry 
a  prior  estimate  of  the  nontidal  sea  surface  heights  [Ray  and 
Byrne ,  2010].  We  have  here  used  weekly  gridded  sea  level 
anomalies  derived  from  a  multi  satellite  analysis  [Pascual 
et  al ,  2006].  The  resulting  tidal  standard  errors,  shown  in 
Figure  lb,  are  mueh  reduced  with  this  technique.  However, 
these  reduced  amplitudes  can  still  reach  1  cm  in  high- 
mesoseale  regions,  whieh  ean  exeeed  internal  tide  ampli¬ 
tudes,  as  is  evident  below. 

3.  Model-Data  Comparisons 
3.1.  Barotropie  Tide 

[14]  Since  the  internal  tides  in  HYCOM  are  generated 
by  the  interaction  of  the  barotropie  tide  with  the  bottom 
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Figure  4.  (a)  M2  surface  tidal  elevation  error  for  the  HYCOM  surface  tidal  elevation  measured  against 
TPXO.  The  contributions  to  the  surface  tidal  elevation  error  resulting  from  errors  in  (b)  tidal  amplitude 
only  and  (c)  amplitude-weighted  phase  following  the  derivation  in  (3).  Units  are  in  centimeters. 
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Figure  5.  Same  as  Figure  4,  but  for  the  K!  tidal  constituent. 
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Figure  6.  M2  internal  tide  amplitude  along  ascending  tracks  from  the  HYCOM  (red)  and  altimeter-based 
analysis  (black).  For  each  track,  the  line  showing  the  coordinates  of  the  track  represents  a  zero  amplitude 
for  the  tides  on  that  track.  The  short-scale  smoothness  is  due  in  part  to  the  application  of  the  band-pass 
filter  and  is  not  due  to  the  response  method  used  in  the  altimetric -based  analysis. 


topography,  an  accurate  barotropic  tide  is  needed  in  order  to 
produce  an  accurate  internal  tide.  To  assess  the  accuracy  of 
the  simulated  barotropic  tide  in  HYCOM,  each  of  its  eight 
tidal  constituents  (computed  from  total  SSH,  which  is  dom¬ 
inated  by  the  barotropic  tide)  are  compared  to  those  from  an 
altimetry  constrained  barotropic  tide  model  (TPX07.2;  an 
update  to  that  described  by  Egbert  et  al  [1994]).  A  recent 
assessment  of  altimeter-constrained  models  [Ray  etaly  201 1  ] 
suggests  M2  RMS  errors  of  about  1 .5  cm  or  less  in  the  deep 
ocean  and  anywhere  from  two  to  ten  times  larger  errors  in 
shallow  water,  depending  on  location.  The  TPX07.2  model 
has  comparable  statistics,  while  nonassimilativc  global  tide 
models  have  much  larger  RMS  errors. 

[15]  Results  from  this  comparison  for  M2  and  Kj  (the 
largest  amplitude  semidiurnal  and  diurnal  constituents)  are 
shown  in  Figures  2  and  3,  respectively.  Qualitatively  the 
tidal  amplitudes  and  phases  in  HYCOM  arc  similar  to  the 


results  from  TPXO,  but  there  are  differences.  One  difference 
is  that  HYCOM  includes  internal  waves,  resulting  in  small 
amplitude,  small  horizontal  scale  perturbations  to  both  the 
amplitudes  and  phases  in  Figures  2b  and  3b.  Another  differ¬ 
ence  is  that  our  HYCOM  tide  simulation  is  a  forward  (non- 
assimilative)  calculation  and  our  barotropic  tides  therefore 
are  not  as  accurate  as  those  in  barotropic  data-assimilativc 
global  tidal  models  such  as  TPXO,  or  in  regional  models 
forced  by  data-assimilative  barotropic  models  at  their  bound¬ 
aries  [e.g.,  Cummins  et  aly  2001;  Merrifield  et  aly  2001]. 

[16]  To  quantify  the  differences  between  HYCOM  and 
TPXO  surface  tidal  elevations,  we  calculate  the  mean  square 
error  (MSE), 

MSE  =  j  \AmCoM^HWOU  -  Atpxo |2,  (1) 
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Figure  1.  The  M2  internal  tide  amplitude  from  the  (a)  alti metric-based  and  (b)  HYCOM  tidal  analyses. 
The  five  subregions  denoted  by  black  boxes  in  Figure  7b  are  used  to  compute  the  area-averaged  ampli¬ 
tudes  in  Table  2. 


where  A  and  <p  are  tidal  amplitude  and  phase,  respectively. 
The  MSE  for  the  eight  constituents  forced  in  the  model  arc 
given  in  Table  1  with  total  MSE  of  9.52  cm  and  MSE  for  the 
leading  semidiurnal  constituent  (M2)  of  7.48  cm  and  leading 
diurnal  constituent  (Ki)  of  2.25  cm.  The  geographical 


distribution  of  the  MSE  for  M2  and  Ki  are  shown  in 
Figures  4a  and  5a.  The  semidiurnal  errors  are  largest  around 
the  continental  margins  and  Southern  Ocean  where  differ¬ 
ences  in  the  bathymetry  of  TPXO  and  our  model  arc  the 
largest  Two  large  regions  of  error  arc  found  in  the  central 
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Table  2.  Area-Averaged  Amplitudes  of  Semidiurnal  Internal 
Tides  From  the  Altimctric-Bascd  (Upper  Value)  and  HYCOM 
(Lower  Value)  Tidal  Analyses  Computed  Over  the  Five  Subregions 
Depicted  in  Figure  7b* 


M2 

s2 

k2 

n2 

Hawaii 

Altimeter 

0805 

0  414 

0  127 

0.209 

HYCOM 

0.887 

0.509 

0  146 

0.162 

Last  of  Philippines 

Altimeter 

0.837 

0.399 

0.121 

0.226 

HYCOM 

0.810 

0.436 

0,160 

0.174 

TropicaJ  South  Pacific 

Altimeter 

0.843 

0.271 

0.083 

0200 

HYCOM 

0.806 

0.283 

0.089 

0.163 

Tropical  SW  Pacific 

Altimeter 

0.758 

0.386 

0.115 

0207 

HYCOM 

0.617 

0.315 

0.111 

0,141 

Madagascar 

Altimeter 

0.715 

0.407 

0  122 

0.192 

HYCOM 

0.665 

0.357 

0.124 

0.128 

Rest  of  world  ocean 

Altimeter 

0.024 

0.242 

0.082 

0.132 

HYCOM 

0,006 

0,001 

0  005 

0.004 

“The  five  subregions  are  hoi  spots  for  generation  of  semidiurnal  tides. 
Unils  are  in  centimeters. 

‘TVenotes  the  area-averaged  amplitude  for  the  rest  of  the  world  ocean 
outside  the  five  hot  spot  regions. 


North  and  South  Pacific.  For  the  diurnal  tides  the  errors  arc 
largest  around  the  continental  margins  and  over  much  of  the 
Atlantic,  Indian  and  Southern  Oceans  The  errors  in  the  tidal 
elevations  can  arise  from  a  combination  of  errors  in  the 
amplitude  A  and  the  phase  <t>.  The  MSE  in  ( 1 )  can  be  rewritten 
as 


MSE  =  -  ((A HYCOM  COS <$> HYCOM  ~  AjPXO  COS 4>tpxo)2 
A-  {A HYCOM  sm(t>HYCOM  “  AfPXO  SUUpTPXo)2^  i 


(2) 


where  MSE  consists  of  contributions  resulting  from  differ¬ 
ences  from  the  in  phase  (cosine)  and  quadrature  (sine)  terms 
for  the  constituent  of  interest.  For  all  constituents,  the  global 
MSE  is  approximately  equally  divided  between  the  in  phase 
and  quadrature  terms.  For  the  semidiurnal  tide,  the  large 
central  North  Pacific  error  is  predominately  in  phase  and  the 
central  South  Pacific  is  predominately  in  quadrature  (maps 
not  shown).  However,  significant  quadrature  errors  in  the 
North  Pacific  and  in  phase  errors  in  the  South  Pacific  are 
found  in  the  same  regions. 

[n]  We  arc  not  just  interested  in  the  generation  ofbaro- 
tropie  and  internal  tides,  but  how  errors  in  the  barotropic  tide 
translate  into  errors  in  the  internal  tides.  As  an  alternative  to 
(2),  we  can  partition  the  MSE  into  contributions  from  dif¬ 
ferences  in  the  amplitude  only  and  from  the  eosinc  of  the 
differences  in  the  phases  weighted  by  the  geometric  mean  of 
the  amplitudes  (amplitude-weighted  phase  errors), 

MSE  =  j j(Ahycom  -  Atpxo)2^ 

A-  [AhYCOmAtpxoO  “  C°S(<P HYCOM  ~  <t>TPXo))]i  P) 


with  the  first  term  on  the  right  hand  side  denoting  a  contri¬ 
bution  to  the  surfaee  tidal  elevation  error  resulting  from 
errors  in  tidal  amplitude  only  (MSE<,mp//fwfc)  and  the  second 
term  (MSE^ltf)  from  errors  in  the  amplitude-weighted 
phase.  As  can  be  seen  from  (3),  if  cither  model  has  small 
amplitude,  then  the  amplitude- weighted  phase  error  will  be 
small  regardless  the  difference  in  phase.  To  illustrate  the 
value  of  this  partitioning,  consider  the  case  of  two  sine 
waves  which  differ  only  in  phase.  From  (2),  the  relative 
contributions  to  the  in  phase  and  quadrature  errors  will  vary 
depending  upon  the  phase  difference.  However,  from  (3), 
the  amplitude  error,  MSE  amplitude  will  be  zero  regardless  of 
the  phase  difference  and  the  error  will  be  only  in  the 
amplitude-weighted  phase  MSE^^*.  Maps  of  the  total 
(MSE),  amplitude  (MSE^//^)  and  amplitude-weighted 
phase  (MSE «w)  errors  for  M2  and  Kj  are  shown  in 
Figures  4  and  5.  The  globally  averaged  statistics  for  all  eight 
constituents  in  HYCOM  arc  listed  in  Table  1  with  the  errors 
approximately  evenly  split  between  amplitude  and  phase. 
However,  unlike  the  in  phase  and  quadrature  errors  from  (2), 
for  the  M2  tide,  the  large  error  regions  in  the  Pacific  are 
predominately  amplitude-weighted  phase  errors.  The  baro- 
tropie  tide  generates  the  baroelinie  tides  through  topographic 
interactions.  These  generation  regions  are  not  uniformly 
distributed  around  the  globe.  In  particular,  large  barotropic 
phase  errors  in  the  Pacific  generation  regions  will  lead  to 
large  baroelinie  phase  errors  associated  with  the  timing  of 
the  generation  of  the  internal  tide.  Thus,  the  MSE  for  the 
internal  tides  will  be  large  due  to  the  timing  errors,  while  the 
amplitudes  will  compare  well. 

[18]  When  the  model  barotropic  tide  is  compared  to  the 
102  pelagic  gauges  described  by  Shum  et  al.  [1997],  the 
RMS  errors  increase  slightly  to  7.80  cm  for  M2  and 
10.22  cm  for  all  eight  constituents,  but  are  still  lower  than  in 
the  work  by  Arbic  et  al  [2004,  2010].  Data  assimilation  for 
TPXO  reduces  the  errors  in  the  shallow  water  tidal  models 
relative  to  the  102  pelagic  gauges  to  ~1.6  cm  for  M2  and 
~3  cm  for  the  eight  constituents  [Shum  et  al.  1997].  When 
compared  to  a  data-assimilativc  model  (TPX07.2),  our 
model  tides  arc  comparable  to  other  nonassimilative  shallow 
water  tide  models,  ~7  em  [Jayne  and  St.  Laurent ,  2001; 
Arbic  et  al,  2004]  and  ^5  cm  for  M2  [Egbert  et  al. ,  2004]. 
Note  that  Jayne  and  St.  Laurent  [2001],  Arbic  et  al.  [2004], 
and  Egbert  et  al.  [2004]  all  utilized  a  rigorous  (i.e.,  non- 
sealar)  SAL  eorreetion,  in  contrast  to  the  HYCOM  results 
shown  here.  Egbert  et  al  [2004]  find  that  5-10%  random 
errors  in  the  bathymetry  can  lead  to  ~8  em  RMS  differences 
in  the  M2  amplitude. 

3.2.  Internal  Tide 

[19]  The  barotropic  tides  interact  with  topography  to 
generate  internal  tides.  Thus,  errors  in  the  barotropic  tides  or 
bathymetry  will  lead  to  errors  in  the  internal  tides.  Globally, 
the  barotropic  tide  errors  are  split  almost  evenly  between 
amplitude  and  phase  errors.  Sinee  phase  errors  in  the  baro- 
tropie  tide  will  cause  phase  errors  in  the  baroelinie  tides, 
the  traditional  RMS  error  statistic  for  the  baroelinie  tidal 
heights  may  not  be  a  good  measure  of  the  model  perfor¬ 
mance.  For  example,  consider  the  M2  internal  tides  in  the 
northeastern  Paeifie,  shown  in  Figure  6,  where  the  amplitude 
of  the  M2  internal  tides  from  the  model  (red)  and  altimeter 
(black)  arc  plotted.  Qualitatively,  the  amplitudes  of  the 
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Figure  8.  The  Kj  internal  tide  amplitude  from  the  (a)  altimctric-bascd  and  (b)  HYCOM  tidal  analyses. 
Areas  where  mcsoscalc  variability  contaminates  the  altimctric-bascd  tidal  analysis  arc  identified  by  the 
red  circles  in  Figure  8a.  The  three  subregions  denoted  by  black  boxes  in  Figure  8b  are  used  to  compute 
the  area-averaged  amplitudes  in  Table  3. 
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Table  3.  As  in  Tabic  2  but  for  the  Three  Subregions  Depicted  in 
Figure  8b,  Which  Arc  Hot  Spots  for  Diurnal  Internal  Tides* 


K, 

0, 

Pi 

Qi 

Central  Indian  Ocean 

Alnmeter 

0.275 

0.185 

0.089 

0.112 

HYCOM 

0.264 

0  149 

0.096 

0.039 

Philippines 

Ahimeter 

0.402 

0.320 

0.125 

0.140 

HYCOM 

0  395 

0.341 

0.153 

0.073 

Central  Tropical  Pacific 

Ahimeter 

0.251 

0.183 

0.081 

0.107 

HYCOM 

0.243 

0  181 

0.117 

0.046 

Resi  of  world  ocean  equaiorward 
of  30° 

Aliimeicr 

0.180 

0.148 

0.060 

0.101 

HYCOM 

0  123 

0.092 

0.056 

0.029 

*Noie  thai  the  average  amplitudes  for  the  diurnal  lides  are  calculated 
equaiorward  of  30°  only.  Uniis  are  in  ceniimeicrs. 


observed  altimetnc  tide  and  model  tide  agree  well  with  RMS 
amplitudes  of  0.361  cm  and  0.346  cm  for  the  altimeter  and 
model,  respectively.  However,  RMS  error  of  the  complex 
amplitudes,  including  the  phase  of  the  tides,  is  0.232  cm.  If 
we  partition  this  difference  into  amplitude  and  phase  errors 
following  (3),  then  the  amplitude  error  is  0. 1 28  cm  while  the 
phase  error  is  0.193  cm.  Thus,  most  of  the  differences 
between  the  model  and  altimeter  internal  tides  arise  from 
phase  errors.  Given  the  errors  in  the  phase  of  the  barotropic 
tide  and  model  bathymetry  errors,  it  is  not  very  surprising 
that  phase  errors  may  dominate  the  internal  tides.  Making 
sure  we  convert  the  proper  amount  of  energy  from  the  bar¬ 
otropic  tide  into  the  baroclinic  tide  is  an  important  first  step 
in  the  evaluation  of  the  model  tides.  We  will  therefore  use 
the  area-averaged  absolute  value  of  the  amplitude  as  the 
statistic  for  our  comparisons.  Using  the  absolute  value  of  the 
amplitude  and  areal  averaging  reduces  the  sensitivity  of 
the  statistics  to  phase  errors. 

[20]  The  global  M2  along-track  altimctric  tidal  analysis 
(Figure  7a)  exhibits  several  internal  tide  generation  regions 
(“hot  spots”)  near  Madagascar,  Hawaii,  east  of  the  Philippines 
and  the  tropical  south  and  southwest  Pacific.  Internal  tides 
radiating  over  long  distances  arc  also  evident,  for  example 
between  the  Aleutian  Islands  and  the  Hawaii  hot  spot  [e.g., 
Cummins  et  a 2001].  Amplitudes  fall  sharply  and  arc  rela¬ 
tively  low  outside  these  hot  spot  regions,  although  close 
analysis  can  reveal  internal  tide  signals  even  in  “quiet”  regions 
such  as  the  southeast  Pacific.  HYCOM  exhibits  similar  fea¬ 
tures  to  those  noted  in  the  altimctric  tidal  analysis  (Figure  7b). 

[21]  To  quantitatively  assess  how  well  the  internal  tide 
results  from  HYCOM  compare  with  the  altimctric-based 
analysis,  arca-avcragcd  amplitude  is  computed  over  five 
subregions  centered  on  internal  tide  generation  regions 
(black  boxes  in  Figure  7b).  In  addition,  area-averaged  sta¬ 
tistics  arc  also  computed  over  the  world  ocean  outside  of 
these  five  hot  spot  regions.  All  four  semidiurnal  constituents 
largely  share  these  hot  spot  regions,  and  summary  statistics 
for  these  constituents  arc  shown  in  Table  2. 

[22]  The  arca-avcragcd  amplitude  is  found  to  agree  well 
across  the  five  hot  spot  subregions  for  the  four  semidiurnal 
constituents  The  average  percent  discrepancy  (() hycom  — 
altim\(\altim\)  x  1 00)  across  the  five  hot  spot  subregions  for  all 
four  constituents  is  ~  1 5%,  with  M2  having  the  lowest  average 
percent  discrepancy  (~9%)  and  N2  having  the  highest 


(~26%).  Across  the  four  semidiurnal  constituents  the  largest 
discrepancy  is  noted  for  the  world  ocean  outside  of  the  five  hot 
spot  regions,  where  the  average  percent  discrepancy  is  ~9I% 
with  the  model  underestimating  the  internal  tide  energy  com¬ 
pared  to  the  altimeter.  Inaccuracies  in  the  simulated  barotropic 
tide,  which  generates  the  internal  tide,  account  for  part  of  the 
discrepancy.  For  example,  the  model  internal  tides  are  too 
weak  in  the  North  Atlantic,  where  the  model  barotropic  tide  is 
weaker  than  the  data-assimilative  barotropic  tide  (Figure  2). 
Another  source  of  the  discrepancies,  mcsoscale  leakage,  will 
be  discussed  later  in  this  section. 

[23]  The  global  Kj  internal  tide  amplitudes  from  the  alti¬ 
mctric  analysis  and  HYCOM  arc  shown  in  Figure  8.  The 
altimetnc  and  HYCOM  tidal  analyses  exhibit  three  main  hot 
spot  regions:  near  the  Philippines,  the  central  Indian  Ocean 
and  the  central  tropical  Pacific.  The  average  percent  dis¬ 
crepancy  across  the  three  hot  spot  subregions  for  all  four 
diurnal  constituents  (Table  3)  is  ~23%,  with  Ki  having  the 
lowest  average  percent  discrepancy  (~3%)  and  Q\  having 
the  highest  (~57%).  Across  the  four  diurnal  constituents  the 
average  percent  discrepancy  for  the  world  ocean  outside  of 
the  three  hot  spot  regions  and  equatorward  of  30°  is  ~37%. 

[24]  Poleward  of  30°  latitude  the  altimctric-based  tidal 
analysis  exhibits  significantly  higher  amplitudes  than  the 
HYCOM  analysis  (Figure  8).  These  high-amplitude  areas 
(circled  regions  in  Figure  8a)  coincide  with  areas  of  high- 
mesosealc  activity,  including  the  Kuroshio,  Gulf  Stream, 
Antarctic  Circumpolar  Current  (ACC)  and  Brazil -Malvinas 
confluence.  This  pattern  is  consistent  across  all  four  of  the 
diurnal  constituents. 

[25]  As  mentioned  in  the  introduction,  propagating  diurnal 
internal  tides  do  not  exist  poleward  of  approximately  30° 
[Gill,  1982].  HYCOM  diurnal  tidal  amplitudes  obtained  from 
hourly  samples  satisfy  this  theoretical  constraint  (Figure  8b). 
However,  as  discussed  in  the  literature  [ Tierney  et  al ,  1998, 
Carrere  et  al ,  2004;  Ray  and  Byrne ,  2010],  the  altimctric 
analysis  (Figure  8a)  shows  features  that  result  from  the  leakage 
of  mcsoscale  activity  into  tidal  frequency  estimates.  This  leak¬ 
age  is  visually  evident  across  all  the  diurnal  constituents,  where 
internal  tides  do  not  propagate,  and  it  can  be  seen  m  the  semi¬ 
diurnal  constituents  as  well.  For  example,  mcsoscale  leakage 
can  be  clearly  seen  in  S2  altimctric  internal  tidal  amplitudes 
(Figure  9a),  with  large  amplitudes  in  the  Kuroshio,  Gulf  Stream 
and  ACC  regions  not  present  in  HYCOM  (Figure  9b).  It  is 
worth  emphasizing,  however,  the  extremely  small  amplitudes 
in  both  Figures  8  and  9.  In  each  case  the  color  bar  spans  only 
5  mm.  It  is  thus  understandable  that  detection  and  mapping  of 
such  small  signals  is  extremely  challenging  for  satellite  altim¬ 
etry,  even  after  almost  two  decades  of  data. 

[26]  Quantitative  evidence  of  mcsoscale  leakage  in  the 
semidiurnal  constituents  is  also  evident  in  Tabic  2,  where 
the  average  percent  discrepancy  over  the  world  ocean  out¬ 
side  of  the  hot  spot  regions  is  80%.  This  discrepancy  is 
significantly  higher  than  the  diurnal  case  (~37%)  because 
the  latter  statistic  was  computed  over  the  30°S-30°N  latitude 
range,  effectively  filtering  out  large  areas  of  mcsoscale 
leakage  (e.g..  Gulf  Stream,  Kuroshio,  ACC). 

4.  Summary  and  Conclusions 

[27]  The  potential  for  the  realistic  simulation  of  barotropic 
and  internal  tides  in  a  high-resolution  global  ocean 
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Figure  9.  The  S2  internal  tide  amplitude  from  the  (a)  altimctnc-based  and  (b)  HYCOM  tidal  analyses. 
Areas  where  mcsoscalc  variability  contaminates  the  altimctric-bascd  tidal  analysis  are  identified  by  the 
red  circles  in  Figure  9a. 


circulation  model  is  examined  using  results  from  year  2006 
of  a  seven  and  a  half  year  1/12.5°  global  simulation  of 
HYCOM  that  resolves  internal  tides,  along  with  barotropic 
tides  and  the  eddying  general  circulation.  Barotropic  tides 
from  HYCOM  are  compared  with  barotropic  tides  from  an 


altimetry-constrained  barotropic  tide  model  (TPX07.2; 
an  update  to  that  described  by  Egbert  et  al.  [1994]).  The 
HYCOM  barotropic  tides  are  comparable  in  amplitude  and 
phase  to  other  nonassimilative  tidal  models.  The  errors  in  the 
HYCOM  barotropic  tide  are  split  almost  evenly  between 
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amplitude  and  phase  errors.  Internal  tide  amplitudes  from 
HYCOM  are  compared  with  results  from  a  TOPEX/POSEI- 
DON  Jason  along-traek  altimetric  tidal  analysis  [Ray  and 
Mitchum ,  1996;  Ray  and  Zaron ,  2011].  This  work  is  the 
first  global  quantitative  comparison  of  the  simulated  internal 
tide  field  to  an  observed  data  set.  The  area-averaged  internal 
tide  amplitudes  over  the  energetie  tidal  regions  are  found  to 
be  quite  similar  to  those  in  along-traek  satellite  altimeter  data 
for  seven  of  the  eight  tidal  constituents  in  HYCOM,  with  Qt 
exhibiting  the  poorest  agreement.  Q]  is  the  smallest  ampli¬ 
tude  ttdal  constituent  and  suffers  from  a  poor  signal-to-noise 
ratio  tn  the  altimeter  data. 

[28]  Away  from  the  hot  spots,  the  comparison  between 
the  model  and  altimctrie  amplitude  is  not  as  good  due,  in 
pan,  to  two  problems,  errors  in  the  model  baro tropic  tides 
and  overestimation  of  the  altimctrie  tides  in  regions  of  strong 
mesoseale  eddy  activity  due  to  leakage  of  this  activity  into 
the  altimetric  tidal  analysis.  This  leakage  affects  all  con¬ 
stituents  and  is  probably  unavoidable,  owing  to  limitations 
tn  time  sampling,  especially  at  the  very  small  (mm  level) 
signal  amplitudes  of  some  internal  tides.  Research  into  fur¬ 
ther  understanding  deficiencies  in  the  barotropie  tides,  and 
hence  the  internal  tides,  in  HYCOM  is  ongoing. 

[29]  These  results  represent  an  encouraging  first  step  in  the 
modeling  of  internal  tides  tn  a  global  ocean  model  that  also 
resolves  the  barotropie  tides  and  eddying  general  circulation. 
This  model,  forced  only  by  atmospheric  forcing  and  the 
astronomical  ttdal  potential,  is  able  to  generate  internal 
waves  over  energetie  tidal  regions  statistically  consistent 
with  observations  without  the  benefit  of  data  assimilation. 
Improvements  to  the  accuracy  of  the  simulated  internal  tide 
arc  likely  to  arise  from  finer  horizontal  resolution,  which 
leads  to  better  resolved  and  represented  oceanic  features  that 
affect  stratification  Reducing  the  errors  in  the  simulated 
barotropie  tide  and  improving  the  accuracy  of  the  bottom 
topography  that  plays  a  key  role  in  internal  tide  generation 
will  also  help.  These  improvements  arc  presently  underway 
in  our  global  HYCOM  development  effort 
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